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First-principles calculations and a tight-binding analysis predict that the iron-pnictide BaCrFeAs2 
is a promising candidate for half-metallic material with fully-compensated magnetization. The 
transition-metal ions Cr and Fe prefer the three-dimensional intervening lattice, which yields the 
antiferromagnetic order of spin orientations. Due to the difference between Cr and Fe in the elec- 
tronegativity, a band gap is opened at the Fermi level in the spin channel in which Fe provides the 
majority carriers. The selective hybridization between 3d orbitals of Cr and As:4p states due to the 
peculiar lattice structure of the iron-pnictide is shown to be crucial for the novel properties. 



I. INTRODUCTION 

As one of the elements well explored from the ancient 
time, iron (Fe) had not been a central player in the mod- 
ern electronic technologies. However, the recent discov- 
ery of superconductivity in iron-pnictide materials pLj cer- 
tainly has opened a new chapter of study on Fe. The 
parent iron-pnictide materials are metallic and antifor- 
romagentic (AFM), and superconductivity is realized by 
carrier doping [1] or application of pressure [2] . So far the 
main focus on the iron-pnictide materials has been su- 
perconductivity, and the AFM order is to be suppressed 
as a competitor. However, it has been revealed during 
the investigation of the mechanism of superconductivity 
that the AFM order of Fe ions itself is also of interest 
which is established by the peculiar lattice structure of 
iron-pnictides via the ferro-orbital ordering [3 . It will 
be fantastic if one can successfully explore the possibility 
of iron-pnictides as a platform for novel spin-dependent 
transport properties in which the magnetic order plays a 
direct role. 

Of many possibilities, we focus here on a novel 
class of materials called half-metallic antiferromagnet 
(HMAFM), which is metallic in one spin channel and 
exhibits a gap at the Fermi level, thus insulating, in 
the other spin channel while showing zero total magne- 
tization in a unit cell [4 . HMAFM can generate spin- 
polarized current without perturbing the surrounding el- 
ements magnetically in a device, and thus is particularly 
useful in spintronics applications. HMAFM can also be 
the parent material for single spin superconductivity with 
triplet Cooper pairs due to the finite density of states 
in only one spin channel [5 . Since the first proposal 

, the possibility of HMAFM has been investigated pri- 
marily in Heulser alloys [6, 7 and double perovskites [8l- 
[T4] . The transition- metal (TM) chalcogenide superlattice 
[15], thiospinels [16] and some disordered systems, such 
as diluted magnetic semiconductors [171 HH]? vacancy- 
induced TM oxides [T^ and hole-doped cuprates [20 are 
also predicted to be the possible condicates. Up to date, 
it is still lacking of a definite confirmation of such novel 
material. 

It is intriguing to notice that the iron-pnictide mate- 



rials can be good candidates for HMAFM. The ground 
state of these materials is AFM and poorly metallic. In 
order to achieve the HMAFM, one only needs to open 
a gap in one spin channel while keeping the AFM or- 
der with compensated magnetization. In the simplest 
picture, a Fe atom has six 3d electrons and shows an 
effective spin moment of 4 /iB due to the intra-atomic 
Hund's coupling. In order to achieve HMAFM, one can 
choose the Chromium (Cr) atom to substitute half of 
the Fe atoms since Cr possesses four 3d electrons. It 
is expected that the Cr atom will not change the AFM 
order of compensated magnetization of the parent mate- 
rial, while modify the band structure due to the different 
atomic number from Fe. 

In the present work, we focus on a typical iron-pnictide 
BaFe2As2 [21 , which belongs to the so-called '122' fam- 
ily. Among the well-studied iron-pnictides, the '122' fam- 
ily has the shortest As- As bond along the c axis, which 
enhances the interlayer coupling. From first-principles 
calculations on the electronic structure of BaCrFeAs2, 
it is revealed that (i) Cr and Fe ions prefer a three- 
dimensional (3D) intervening lattice; (ii) the magnetic 




FIG. 1. (a) Crystal structure of the TM-As layer of 
BaCrFeAs2; (b) Spin-density of the AFM state (blue/red sur- 
face: spin-up/spin-down). Isovalue=0.025 e/A^. 
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moments of Cr and Fe ions are aligned antiferromagnet- 
ically, and with small but non-negligible contributions 
from As atoms parallel to that from Fe ions, the total spin 
magnetization in a unit cell is zero; (iii) a gap is opened 
at the Fermi level (£^f) in the spin channel where Fe elec- 
trons contribute as majority carriers by the hybridization 
between As:4p states and minority Cr:3d states due to 
the difference in the electronegativity between Fe and Cr 
as well as the peculiar crystal structure of iron-pnictide. 
Based on these results, we predict BaCrFeAs2 as a pos- 
sible HMAFM. 

This paper is organized as follows. In the next section 
the computational details of calculations are briefly de- 
scribed. In the first part of section III, we focus on the 
electronic structure of BaCrFeAs2; then a tight-binding 
(TB) picture is presented. Finally the summary of the 
paper is given in section IV. 



II. COMPUTATIONAL DETAILS 

First-principles calculations were performed by VASP 
package [22j [23 within the scheme of PAW method 
[24] and generalized gradient approximation [25^ for 
the exchange-correlation functional [26]. 8x8x4 
Monckhorst-Pack special points and 500 eV are used for 
k sampling and plane- wave basis set. Starting from the 
experimental lattice parameters of BaFe2As2 in tetrago- 
nal structure [27 , both the ion coordinations and lattice 
parameters are fully relaxed until the force and stress are 
less than 0.01 eV/A and 1 kBar. Calculations were also 
performed by the Quantum-Espresso package [28 , and a 
good agreement is achieved. 

The on-site energies and hopping integrals of the TB 
model are calculated by the maximally localized Wannier 
functions (MLWF) method as implemented in the Wan- 
nier90 package [29] jointly with Quantum-Espresso [28] . 
The Bloch wavefunction of BaCrFeAs2 is first calculated 
in density functional formalism. By specifying the 3d or- 
bitals of TM ions and 4p orbitals of As ions as the initial 
guess, the periodic wavefunctions are transformed to the 
MLWF representation in real space. Meanwhile, the on- 
site energies and hopping integrals can be obtained. It 
is noticed that the TB model including only the nearest 
neighboring hoppings can reproduce well the DFT band 
structures. 



III. RESULTS AND DISCUSSION 

A. First-principles results 

In BaFe2As2, Fe ions tetrahedrally coordinated by As 
atoms form Fe-As layers, which are sandwiched by Ba 
layers. We concentrate on the situation that half of the 
Fe atoms are substituted by Cr atoms in a x ^/2 x 1 
cell. The possible atomic configurations which have been 
studied include: (1) cross-checkerboard type as shown in 
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FIG. 2. (a) Total density of states (DOS) of half-metallic 
antiferromagnet BaCrFeAs2. (b) and (c) Projected DOS of 3d 
states for Cr and Fe, respectively. The Fermi energy is set 
to energy zero. The bonding, anti-bonding and non-bonding 



states of t2g are labeled as t 
Projected DOS of As:4p states. 



t2g and t§^. 



respectively, (d) 



Figure lb, (2) checkerboard type with same positions of 
Cr and Fe ions along c axis, (3) alternating Cr and Fe 
layers, and (4) alternating Cr and Fe stripes along a or 
b axis. We calculate the total energies of them with all 
the possible spin orderings of the TM ions. It is revealed 
that the ground state of BaCrFeAs2 is characterized by 
the intervening Fe and Cr ions and 3D AFM spin order, 
as shown in Figure lb. The resultant lattice parameters 
of a=5.725 A and c=13.290 A are in good agreement with 
the experiments [30 . 

The magnetization is mainly contributed by Fe (-2.62 
/ie) and Cr (2.75 /ie) ions. In addition, each As atom car- 
ries a small magnetic moment (-0.08 jub) parallel to that 
of Fe ions due to the hybridization between As:4p and 
TM:3d electrons. The Ba layers act as charge reservoir, 
and thus the magnetic moment of Ba^+ ions is negligible. 
According to the calculation, the magnetic moments of 
the Fe, Cr and As atoms compensate completely, result- 
ing in the total magnetization of 0.0 fiB per unit cell. 

The total density of states (DOS) of BaCrFeAs2 is de- 
picted in Figure 2a. While the spin-up channel exhibits 
an evident metallic character similar with the parent ma- 
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terial [31], a gap of ~ 0.3 eV emerges at in the spin- 
down channel. Combining the zero magnetization and 
half-metahic band structure, BaCrFeAs2 is expected to 
be a HMAFM. 

Let us look into the detailed electronic band structure 
of the new material BaCrFeAs2. In the spin- up channel, 
Cr nominally has four 3d electrons, and thus its majority 
3d band is not fully occupied, as illustrated by Figure 2b. 
The hybridization between As:4p and Cr:t2^ states forms 
the bonding states (t2^) deep in the valence band, and 
partially occupied anti-bonding states (tf^) around E^. 
Although there exist Fe:3d-As:4p as well as Fe:3da;2_^2- 
Cr:3da.2_^2 hybridizations, most of the Fe:3d minority 
states remain non-bonding and locate around {t2g and 
d3^2_^2 in Figure 2c). BaCrFeAs2 is therefore metallic in 
the spin-up channel. 

In the spin-down channel, Fe:3d majority states lie 
deeply in the valence band and are fully occupied, as 
shown in Figure 2c. Most of the minority states of Cr 
are above Ep and form the bottom of conduction band, 
in contrast with Fe in the spin-up channel, due to the 
difference in the electronegativity between Cr and Fe. 
Around the top of valence band there are noticeable con- 
tributions from Cr:3d minority states as indicated in Fig- 
ure 2b. It is intriguing to notice that the lattice struc- 
ture of the present iron-pnictide material is crucial for 
opening a gap of ~ 0.3 eV between Cr:3d states in the 
spin-down channel. Since the Ci'.Sdxy orbitals reside on 
the TM plane with their petals pointing to the nearest 
neighboring Cr ions along (110) or (1-10) direction, and 
the As atoms just locate above the center of the TM 
squares (see inset of Figure 3a), the As'Apz orbitals medi- 
ate the exchange interaction between the Ci'.Sdxy states. 
Meanwhile, the As'Apy (As:4pa;) state provides an evi- 
dent bridge for the exchange coupling between the two 
Ci'.Sdxz {Ci'.Sdyz) states. As the result, the spin-down 
As:4p states in the vicinity of E^ hybridize with a part 




TABLE L Hopping integrals between the orbitals of nearest- 
neighboring As- As, As-Cr/Fe and Cr-Fe ions derived from the 
MLWF calculation in the units of eV 



orbitals 



hopping integrals 
Spin-up Spin-down 



As:px /py-Ci-.dyz /dxz 


0.70 


0.80 


As:p2-Cr:dcc^ 


0.61 


0.68 


As:px/Py-Cr:d^2_y2 


0.50 


0.65 


As:p2-Cr:d322_^2 


0.50 


0.65 


As:pcc/py-Cr/Fe:d3^2_^2 


0.24 


0.27 


As:px /Py-Fe:dyz / dxz 


0.66 


0.58 


As:pz-Fe:dxy 


0.57 


0.50 


As:px/py-Fe:d^2_y2 


0.50 


0.39 


As:pcc/py-As:pcc/py (I) 


0.27 


0.27 


As:pa:-As:py (I) 


0.34 


0.33 


As:p^-As:p^ (II) 


0.29 


0.28 


As:pz-As:px/Py (II) 


0.41 


0.38 


As:px/py-As:px/py (II) 


0.34 


0.30 


As:p;2-As:p;, (III) 


0.73 


0.76 


Cr/Fe:da,2_^2-Fe/Cr:d322_^2 


0.22 


0.21 


Ci-.dxz/dyz-Fe-.dxz/dyz 


0.28 


0.28 


Cr:d^2_^2-Fe:d^2_^2 


0.37 


0.37 


Cr-Axy-Fe-Axy 


0.27 


0.26 



^ The Roman numbers denote the hoppings between As- As 
ions (I) within the same As plane, (II) crossing the TM plane 
and (III) along the c axis, which can be read from Figure 1. 



TABLE 11. Onsite energy of As:4p and TM:3d orbitals for the 
TB model derived from the MLWF calculation in the units of 
eV 

Cr:d3^2_^2 CrAxz/dyz Cr:d^: 



Spin up -1.61 
Spin down 0.18 
Fe:d3,2_ 



-1.08 
0.61 



-1.62 
0.20 



,2 FeAxz/dyz Fe:d^2_^2 



Cl Axy 

-1.38 
0.65 

FQAxy 



Spin up -0.24 -0.17 -0.46 -0.14 

Spin down -2.66 -1.96 -2.71 -2.00 

As:p^ As:px/py 



Spin up -2.24 -2.39 
Spin down -2.23 -2.34 



of the Cr:t2^ minority states (Figure 2b and 2d), and 
lower their energies to the valence band. These features 
are well captured by the band-decomposed charge den- 
sity [22j |23] depicted in Figure 3 obtained by calculations 
based on the a/2 x a/2 x 1 cell. 



FIG. 3. Band-decomposed charge density for the (a) Cv.Sdxy 
bonding via As:4p2, (b) Cr:3dx2 bonding via As:4py and (c) 
Cr-.Sdyz bonding via As:4pcc at the F point of the spin-down 
valence band maximum associated with a a/2 x a/2 x 1 cell. 
Only one TM-As layer is shown. Isovalue= 0.004 e/A^ for (a) 
and 0.008 e/A^ for (b) and (c). Insets schematically show the 
corresponding atomic orbitals of Cr and As. 



B. Tight-binding Model 

For a better understanding on the novel HMAFM 
property of BaCrFeAs2, in which many orbitals con- 
tribute around E^^ it is helpful to build a TB picture. For 
this purpose, we adopt the MLWF method ^281 [29], which 
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FIG. 4. Band structures of BaCrFeAs2 with the red sohd and green dotted curves denote results based on the first-principles 
calculations and the TB model, respectively, (a) /(c) for spin- up/down channel; (b) for on-site energies of relevant orbitals of 
the TB model derived from the MLWF calculation. 



yields both on-site energies and hopping integrals of the 
TB model. Since the whole TM:3d and As:4p orbitals 
are involved in the covalency, all of them are treated as 
the initial guess for orbital projections. It is confirmed 
that the DFT results can be perfectly reproduced. We 
then reduce the hopping ranges to the minimum with 
which one can still capture the main features of the DFT 
results. In this way, the important orbitals which govern 
the band dispersion around are identified as TM:t2^ 
and As:4p, and the hopping integrals are for the nearest 
neighbors only. 

The numerical results for the hopping integrals and the 
on-site energies are given in Table I and II respectively, 
and the band structures derived from the TB model are 
shown in Figure 4 in fairly good agreement with the 
VASP results. The onsite energies of spin-up Fe:3d states 
are slightly below E-p^ while the spin-down Cr:3d orbitals 
are slightly higher than E^. It is clear that this difference 
is responsible for the half metallic feature of the present 
material BaCrFeAs2. 

As listed in Table I, the Asip^-Asip^ hopping along c 
axis is quite large (~ 0.7 eV) due to the crystal struc- 
ture of the present '122' system with the shortest As- 
As bond along the c axis. This interaction is crucial 
for the c axis AFM order of BaCrFeAs2, and thus the 
feature of HMAFM, in contrast to the case of the so- 
called '1111' system Therefore we predict that the 
present BaCrFeAs2 as well as other '122' materials are 
the most promising candidates for HMAFM among the 



iron-pnictides. 

IV. CONCLUSIONS 

To summarize, based on first-principles calculations 
and a TB model we predict that BaCrFeAs2 is a can- 
didate material for half- metallic antiferromagnet. The 
novel half-metallic band structure of this material is gen- 
erated by the selective, strong hybridization between 
As:4p states and minority Cr:3d states, due to the differ- 
ence in the electronegativity between Fe and Cr. Both 
the density-functional calculations and the TB analysis 
reveal clearly that the crystal structure of iron-pnictide 
is crucial for opening a gap around as large as ^ 0.3 
eV in the single spin channel. The zero magnetization is 
achieved by the magnetic compensation of Cr, Fe and As 
atoms. It is highly anticipated that the theoretical pre- 
diction formulated in the present work can be checked 
experimentally. The endeavor will then open a new field 
of iron-based spintronics. 

ACKNOWLEDGMENTS 

The authors thank I. V. Solovyev, B. Liu and Y.-M. 
Nie for discussions. The calculations were performed 
on Numerical Materials Simulator (SGI Altix) of NIMS. 
This work was supported by WPI Initiative on Materials 
Nanoarchitronics, MEXT, Japan. 



[1] Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, [3] C.-C. Lee, W.-G. Yin, and W. Ku, Phys. Rev. Lett. 103, 

J. Am. Chem. Soc. 130, 3296 (2008). 267001 (2009). 

[2] M. S. Torikachvih, S. L. Bud'ko, N. Ni, and P. C. Can- [4] H. van Leuken and R. A. de Groot, Phys. Rev. Lett. 74, 

field, Phys. Rev. Lett. 101, 057006 (2008). 1171 (1995). 



5 



[5] W. E. Pickett, Phys. Rev. Lett. 77, 3185 (1996). 

[6] S. Wurmehl, H. C. Kandpal, G. H. Fecher, and C. Felser, 

J. Phys.: Condens. Matter 18, 6171 (2006). 
[7] I. Galanakis, K. Ozdogan, E. §asioglu, and B. Aktas, 

Phys. Rev. B 75, 172405 (2007). 
[8] W. E. Pickett, Phys. Rev. B 57, 10613 (1998). 
[9] J. H. Park, S. K. Kwon, and B. I. Min, Phys. Rev. B 65, 

174401 (2002). 

[10] M. S. Park and B. I. Min, Phys. Rev. B 71, 052405 

(2005) . 

[11] Y. K. Wang and G. Y. Guo, Phys. Rev. B 73, 064424 

(2006) . 

[12] M. Uehara, M. Yamada, and Y. Kimishima, Sohd State 

Commun. 129, 385 (2004). 
[13] Y. K. Wang, P. H. Lee, and G. Y. Guo, Phys. Rev. B 80, 

224418 (2009). 

[14] V. Pardo and W. E. Pickett, Phys. Rev. B 80, 054415 
(2009). 

[15] M. Nakao, Phys. Rev. B 77, 134414 (2008). 

[16] M. S. Park, S. K. Kwon, and B. L Min, Phys. Rev. B 64, 

100403(R) (2001). 
[17] H. Akai and M. Ogura, Phys. Rev. Lett. 97, 026401 

(2006). 

[18] L. Bergqvist and P. H. Dederichs, J. Phys.: Condens. 
Matter 19, 216220 (2007). 

[19] D. Kodderitzsch, W. Hergert, Z. Szotek, and W. M. Tern- 
merman, Phys. Rev. B 68, 125114 (2003). 

[20] Y.-M. Nie and X. Hu, Phys. Rev. Lett. 100, 117203 
(2008). 

[21] M. Rotter, M. Tegel, D. Johrendt, L Schellenberg, 
W. Hermes, and R. Pottgen, Phys. Rev. B 78, 020503(R) 



(2008) . 

[22] G. Kresse and J. Hafner, Phys. Rev. B 47, 558 (1993). 
[23] G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 
(1996). 

[24] P. E. Blochl, Phys. Rev. B 50, 17953 (1994). 

[25] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

[26] It is widely accepted that the electronic structure of iron- 
pnictides (particularly the Fermi surface) predicted by 
the local density approximation (LDA) or GGA is in good 
agreement with the ARPES experiments. In the current 
study we focus on the half-metallic bands in the vicinity 
of the Fermi level. It is expected that GGA can predict 
the accurate band structures of BaCrFeAs2 and therefore 
the strong correlation effect (Hubbard U) is unnecessary. 

[27] Q. Huang, Y. Qiu, W. Bao, M. A. Green, J. W. Lynn, 
Y. C. Gasparovic, T. Wu, G. Wu, and X. H. Chen, Phys. 
Rev. Lett. 101, 257003 (2008).^ 

[28] http : //www . quantum- espresso . org 

[29] A. A. Mostofi, J. R. Yates, Y.-S. Lee, I. Souza, D. Van- 
derbilt, and N. Marzari, Comput. Phys. Commun. 178, 
685 (2008). 

[30] A. S. Sefat, D. J. Singh, L. H. VanBebber, 
Y. Mozharivskyj, M. A. McGuire, R. Jin, B. C. Sales, 
V. Keppens, and D. Mandrus, Phys. Rev. B 79, 224524 

(2009) . 

[31] D. J. Singh, Phys. Rev. B 78, 094511 (2008). 

[32] M. Nakao, J. Phys.: Conf. Ser. 150, 052182 (2009). 



